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SUMMARY 

Analysis is made of the flow field produced by oblique impingement 
of weak plane disturbances of arbitrary profile on a plane normal shock. 
Three types of disturbance are considered: 

(a) Sound wave propagating in the gas at rest into which the shock 
moves. The sound wave refracts either as a simple isentropic sound wave 
or as an attenuating isentropic pressure wave, depending on the angle 
between the shock and the incident sound wave. A stationary vorticity 
wave of constant pressure appears behind the shock. 

(b) Sound wave overtaking the shock from behind. The sound wave 
reflects as a sound wave-, and a stationary -vorticity wave is produced. 

(c) An incompressible vorticity wave stationary in the gas ahead of 
the shock. The incident wave refracts as a stationary vorticity wave, 
and either a sound wave or attenuating pressure wave is also produced. 

Computations are presented for the first two types of incident wave, 
over the range of incidence angles, for shock Mach numbers of 1, 1.5, 
and oo. 


INTRODUCTION 

The unsteady one -dimensional interaction of normal shock waves and 
disturbances, such as sound waves or other shock waves, has been studied 
quite thoroughly (an example is Kantrowitz' paper on shock stability, 
reference l) . The steady interaction between normal shock waves and 
plane Mach waves has been treated by Adams (reference 2). 

The general class of unsteady flow problems is currently of 
increasing interest, in connection particularly with stability of high- 
speed aerodynamic and combustion processes. The effect of a shock pas- 
sing through a flow field (or vice versa) is likely to be important in 
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many applications. For example, a hot-wire anemometer intended to meas- 
ure the fluctuating field of turbulence in a supersonic stream will 
actually measure the turbulence as modified by passage through the 
nearly steady bow shock of the probe. 

Considering, for simplicity, that the flow interacting with a nor- 
mal shock is a nonviscous field of weak disturbance, it may usually be 
considered irrotational and isentropic (such as produced by a moving 
slender body) and therefore can be imagined to be composed of a suitable 
array of sound waves. Another possible type of weak nonviscous distur- 
bance would be a stationary, incompressible flow of variable vorticity 
(turbulence which is convected rapidly past the point of observation is 
commonly thought of in this way) . 

Either of these two types of flow may be represented as a linear 
combination of plane waves (each wave either a sound wave or a rota- 
tional wave, depending on the type of flow to be represented) of various 
amplitudes, wave lengths, and orientations. Thus, if the interactive 
effect of a shock and each constituent wave may be found by a linear 
analysis, the complete problem may in principle be solved by linear com- 
bination of the resulting flow fields behind the shock. The interaction 
between a turbulent field and a wind-tunnel screen and/or contraction 
has been successfully carried out in references 3 through 5 by this 
method . 

The present report concerns the interaction of a normal shock met 
obliquely by a plane sound wave, or by a convected plane vorticity wave. 
Since sound waves may impinge on a shock either from upstream or down- 
stream, both cases are considered. The oblique interaction of a shock 
and weak vorticity wave is also treated in a current investigation by 
Ribner (reference 6). 

The shock is considered to be moving freely into gas nominally at 
rest (as in a shock tube, when wall effects are neglected) . Of course, 
if the observer moves at a constant speed with the shock the flow 
appears as that associated with a steady shock, under different stagna- 
tion conditions. The shock- tube point of view is adopted in order that 
there be no question of how the equivalent steady shock is "anchored"; 
that is, end effects on the shock are not contemplated. 


GOVERNING UNSTEADY EQUATIONS 


In the following paragraphs, the 
pertain to the propagation of a plane 
rest, as modified by the influence of 
disturbance. 


equations will be derived which 
normal shock wave through a gas at 
a weak pattern of unsteady 
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If the shock propagates without disturbance, its instantaneous 
position is x^_ = Vt, in a coordinate system fixed in th'e fluid nomi- 
nally at rest ahead of the shock (fig. l). (All symbols are defined in 
appendix A.) The constant velocity of the shock front is V, and the 
corresponding constant velocity of the gas behind the shock is U, in 
this system. This one-dimensional motion is considered to be perturbed 
slightly by the presence of a weak field of unsteady plane flow. The 
velocity ahead of the shock is written as u 1 (x,y,t), v 1 (x,y,t)} behind 

the shock as U + u 2 (x,y,t), v 2 (x,y,t). Pressure, density, and temper- 
ature are written as P + p, R + p, and 0 +0, respectively, where the 
capitalized symbols refer to the basic steady- shock motion, and the 
lower case to the unsteady disturbance. Throughout, subscripts 1 and 2 
are used to specify conditions ahead of and behind the shock, respec- 
tively. As a result of the unsteady disturbance, the shock front 
itself undergoes a small unsteady displacement, given by £(y,t). Thus, 
the position of the shock at any instant is Vt + £(y,t). 


Shock Relations 

Because of the rapidity with which changes occur across a shock 
wave, the disturbed shock may be regarded as behaving in a locally 
quasi- steady manner j that is, in a coordinate system fixed in the shock 
at each instant, the usual Rankine-Hugoniot relations apply. Because 
of the disturbance, the shock is slightly oblique in such a coordinate 
system (see sketch). 


y 



Because the shock is only slightly oblique, the shock relation concern- 
ing the product of velocity components normal to the shock front, ahead 
of and behind the shock (reference 7) may be written approximately as 
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(V - U + ^ - u 2 ')(V + | t - u x ) = 2 Lzl | (V + | t 



(la) 


The assumption of a slightly oblique shock also provides that the equa- 
tions of conservation of normal energy and mass, respectively, may be 
written: 


| (V + It - u i) 2 + pri (®1 + ©i) = ! (V - u + s t - U 2 ) 2 + (02 + e 2 ) 

(ib) 

(Hi + Pi) (V + | t - up = (R 2 + p 2 )(V - U + | t - u 2 ) (lc) 

The remaining oblique shock relation states that the velocity component 
parallel to the shock is unaltered by passage through the shock. Because 
the shock is assumed to be only slightly oblique, this relation is, 
approximately, 


or, 


V 1 - Vf y = v 2 - (V - U)| y 

\ 


Vo - V-i 

-ty = — — (Id) 

The shock relations for the basic undisturbed shock propagation are 
obtained from equations (la), (lb), and (lc), neglecting small 
quantities : 



(2a) 

2 V 2 * r I J 1 e i - 2 (V - U) 2 + r r J ! 02 

(2b) 

R X V = R 2 (V - U) 

(2c) 


Terms of equations (la), (lb), and (lc) which are of first order in 
small quantities yield the conditions which the disturbance field must 
satisfy at the shock: 


( v " u )(lt “ u l) + v Ut “ u 2) 


= 2 


r - i 
r + i 



Up + 


r J 


r - i 



(3a) 
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v(| t - u i) + e i = < v - u )Ut - u 2) + 7^1: 0 2 

Pl V + RiUt - u l) = P2( V ' u ) + R 2<lt • u 2) 

Equations (2), the state equation , ~ 

P = JR © (basic flow) 

> 

£ = •§ + -^ (perturbation flow) 

r K fc) 


(3b) 

(3c) 


(4) 


and the assumption that the incident flow ahead of the shock is isen- 
tropic permit equations (3) to be simplified, yielding the following set 
of disturbance shock relations (equation (id) is also included), relating 
the disturbance fields ahead of and behind the shock: 


t 


t - u 2 . „ It - U 1 

— — = B 1 v 


^2 _ p It " U 1 
r 2 1 V 



= II Ut 

v v " v 


+ B 





(5a) 

(5b) 

(5c) 

(5d) 


where the coefficients are constants depending on the Mach number of the 
undisturbed shock M = V/a^: 


\ 


/ 


\ 
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b r - i 2 

2 r + 1 rM 2 


C n = 


1 1 + r - ~ 1 M 2 


c 2‘. = 7 




( 6 ) 


Di = 


2rMf 


rM' 


;2 _ r - i 


Do = 




rM' 


2 _ r - i 

2 


J 


Equations of Plane Disturbance Field 


In addition to the perturbation shock relations (equations (5)) 
which are concerned with the compatibility at the shock of weak distur- 
bance fields ahead of and behind the shock, the equations satisfied by 
the disturbance fields themselves are required. The nonviscous equa- 
toins of motion are written in a coordinate system at rest relative to 
the gas ahead of or behind the shock. That is, the following equations 
apply in the coordinate system of figure 1 for the disturbance field 
ahead of the shock, and in a coordinate system moving with velocity U, 
for the disturbance behind the shock. Subject to interpretation of the 
coordinate system, the same . equations apply in both regions, and there- 
fore the subscripts 1 and 2 are omitted for the time being. To first 
order in small quantities: 


Momentum: u^ = - —^p x 

V t - - I Py 

Continuity: + R(u x + v ) 


= 0 


( 7 ) 

\ 

( 8 ) 



\ 


MCA TN 2879 


Energy: C-yRS-f. + P(u x + v y) = 0 


(9) 


State- £ - £ + 5 
State . p - R + 0 

These equations may he combined, to 
disturbance: 

d 2 2 / d 2 

_st 2 " a U 2 

the equation for the entropy disturbance (which may be shown to be pro- 
portional to (p/p) - t(p/p)): 

and the equation for the vorticity (V x - Uy) 

5E < V X - "y> = 0 

Thus, the pressure disturbance satisfies the wave equation (10), and any 
entropy variation (equation (ll)) or vorticity variation (equation (12)) 
is steady, relative to the main flow. Any disturbance field satisfying 
these linear equations may be regarded as composed of two parts, one 
steady, and the other unsteady, in a coordinate system at rest in the 
mai n flow. From equations (ll) and (12), variation of entropy and vor- 
ticity may be assigned to the steady flow and from equation (8)^ the 
associated velocity components satisfy the incompressible continuity 
equation. From equations (7), pressure variations must be assigned to 
the unsteady flow, satisfying the wave equation (10). The unsteady por- 
tion of the flow may then be regarded as produced by a pattern of sound 
waves. A weak nonviscous disturbance field may therefore be considered 
to include: 

1. An unsteady, isentropic, irrotational disturbance, which may be 
regarded as produced by a pattern of sound waves, and 

2. A steady rotational disturbance of constant pressure and (in 
general) variable entropy and density. 

TYPES OF INITIAL DISTURBANCE CONSIDERED 

The present analysis concerns the interaction of a shock wave with 
three types of initial plane disturbances: 


( 11 ) 

of the disturbance flow: 

( 12 ) 


( 4 ) 


yield the equation for the pressure 


£>y 2 J 


p » 0} a 2 = ^ 


( 10 ) 
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A. Sound Wave Overtaken by Shock 

The shock moves into a region in which a plane sound wave is propa- 
gating in a direction oblique to the direction of shock propagation 
(fig. 2(a)). Since the shock velocity is supersonic relative to the gas 
ahead, it will overtake the sound wave, whatever its direction of propa- 
gation. The solution for the interaction of such a wave with a shock 
may in principle be generalized by linear superposition to provide anal- 
ysis of the passage of a shock through any isentropic field of small 
disturbance. 


A general plane sourid wave may be represented as follows (the par- 
ticular profile of the wave need not be specified): 


U 1 

T = A i f 


v i 

T = V 


c 


mxi - Zy + a-|_t' 
*1 

’mxi - Zy + a^t 
*1 


i) 


p, /mx - Zy + a t\ 

Pi = A 3 f ( XI j 

Pi (i 

T ± =^ f \ 


fmx.± - Zy + a^t 




J 


(12a) 


where 


l = sin )|r^} m = cos ^ 


(13) 


and is a length characterizing the scale of the disturbance. If the 

function f were a sine wave, would be equivalent to the wave 
length. By equations (12), (ll), and (7), respectively. 


ZAn = - mA. 

-L Cj 

A 3 = rA 4 


A 


1 


m 

‘ 7m A 3 


"v 

> 


J 


(12b) 
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The disturbance is of the type 1 discussed in the previous section 
(unsteady, isentropic, irrotational) , and is longitudinal; that is, the 
fluctuating velocity component is in the direction of propagation of the 
sound wave. 


B, Shock Overtaken by Sound Wave 

The sound wave propagates relative to the fluid behind the shock, 
in such a mann er as to overtake the shock (fig. 2(b)). Thus, considera- 
tion will be restricted to cases for which -ag cos +21 - U. The 

initial disturbance may be specified in a manner similar to that 
employed for the preceding case. 


The subscripts 1 and 2 have been introduced to denote the flow 
ahead of and behind the shock, respectively. In the present problem, 
the entire flow disturbance occurs behind the shock. The subscript 2 
is therefore appropriate to both the incident and reflected waves, which, 
will hereinafter be distinguished by second subscripts 1 and 2, 
respectively. 


u. 


21 


= An f 


mxg - Zy + 


'21 A _ / mx 2 - 
~ = A 2 1 


y + 

21 ' 


P21 (^z ■ + a 2 ^ 

*T ■ A3f 


^ 21 


p 2 i /mx 2 - Zy + a 2 t 

’ A(lf \ 


''v 


(14a) 


J 


Pertinent equations of motion provide, as before, 

ZA-^ = - mAg ^ 

a 3 = rA 4 


a 2 

A 1 = - m Tv A 3 


(14b) 
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The coordinate system X£,y is fixed relative to the flow behind the 
shock. 


C. Stationary Vorticity Wave Overtaken by Shock 

The shock moves into a region occupied by a stationary plane shear 
disturbance of constant density, oblique relative to the shock front 
(fig. 2(c)). A system of such waves may be employed to represent a tur- 
bulent field (references 3 and 4). Therefore, the effect of the passage 
of a shock through a single oblique shear wave may in principle be gen- 
eralized by Fourier superposition to provide an analysis of the passage 
of turbulence through a normal shock. 

The incident vorticity wave, of arbitrary profile, may be repre- 
sented as follows: 


(15) 



From equation (9), 

mA x = IA 2 (16) 

The disturbance is therefore a special case of the type 2 discussed 
in the previous section (a steady vorticity disturbance of constant pres- 
sure), and, by continuity, must be transverse; that is, the fluctuating 
velocity component is parallel to the plane of the shear waves. Since 
the wave is transverse, there may be a component of velocity disturbance 
parallel to the shock (perpendicular to the plane of fig. 2(b)) which 
may be of arbitrary amplitude . 

This type of interaction (problem C) is treated by Ribner in 
reference 6 . 



ANALYSIS OF INTERACTION BETWEEN SHOCK AND IMPOSED DISTURBANCE 

Three Types of Refraction at Shock 

It has previously been shown that the equations of motion imply 
that any weak disturbance field may be divided into two parts: An 

unsteady isentropic irrotational field; and a steady vorticity 
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disturbance. Thi6 point of view may be adopted with regard to the 
disturbance downstream of the shock, produced by the interaction. 

Simple refracted sound wave - problem A . - In the case of a sound 
wave overtaken by a shock (problem A of the previous section), it would 
seem reasonable to expect that for typ near either 0 or * (fig- 2(a)), 

the isentropic part of the downstream field would simply be a refracted 
sound wave traveling away from the shock. This is indicated by the 
sequence of events shown in figure 3(a). At time tp, the initial wave 

intersects the shock front at point Pp ■ At a later time tg = t-^ + St, 
the sound wave has moved a distance ap&t, the shock has moved a dis- 
tance V5t, and the intersection occurs at point P 2 . In the meantime , 
a cylindrical sound wave has been generated at point Pp as a result of 
the shock interaction and expands with velocity a 2 , while being con- 
vected with a velocity U. Thus, at time t 2 , the effect of the inter- 
section at tp is felt within a cylinder of radius a 2 St, with center 
at point Q 2 . According to figure 3(a), an envelope is formed, and may 
be identified as a simple refracted sound wave. 

Attenuating refracted pressure wave - problem A . - Figure 3(a) is 
drawn for a rather small value of \|fp. If ^p is increased, there 

appears a critical value i|r ci (fig. 3(b)) beyond which no envelope may 
be drawn. Thus, when +p>i C 2> ^e influence of intersection is 

felt at P 2 before the intersection arrives at P 2 . However, as i|fp is 
further increased, there appears another critical angle \|r cu beyond 

which simple envelopes may again be drawn (fig. 3(d)) and simple sound 
wave refraction occurs. 

When <'l r i< ♦cu* downstream pressure disturbance cannot be a 

simple sound wave. The cylindrical sound waves produced by the inter- 
action at the shock do not coalesce, but rather continue to expand 
independently, thus dimini shing in strength as time progresses. Accord- 
ingly, the isentropic part of the downstream disturbance may be expected 
to die out at large distances downstream of the shock. This attenuating 
disturbance may, however, be expected to remain planar, because both the 
incoming disturbance and the shock are plane. This attenuating wave has 
been called a pressure wave rather than a sound wave, because, as will 
be shown subsequently, it does not propagate at the local velocity of 
sound . 
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Expressions for and t cu in terms of shock Mach number may- 

be obtained from the following equation derived by inspection of 
figure 3(b): 


a 2 2 - (V - U) 2 = (V cot \|r c + esc \|r c ) 2 (17) 

The solution of this equation is shown in figure 4, labeled "sound wave". 
The curves labeled "stationary vorticity wave" will be discussed subse- 
quently. The curves both approach a value of 180° for M = 1, and have 
a half -order singularity there. As M->®, the curves become symmetric 
about = 90°, because a^ becomes insignificant compared to V. The 

limiting value of \|r c j -is 67.8°. 

Steady vorticity wave - problem A . - Returning to figure 3, if a 
vorticity disturbance is created at the shock-disturbance intersection, 
and is thence convected with a velocity U, a vorticity wave appears 
along the line connecting P 2 and Q 2 , whatever the value of 

Thus, of the three types of refractions discussed, the "steady vor- 
ticity wave" always appears, in combination with either a "simple 
refracted sound wave", or an "attenuating pressure wave", depending on 
the angle \|r^. 

Expectations for problem B . - In the case of the shock overtaken by 
a sound wave, it may be shown by constructing figures similar to fig- 
ure 3 that a simple reflected sound wave will always occur, in conjunc- 
tion, of course, with a steady vorticity wave. 

Expectations for problem C . - In the case of the shock overtaking a 
stationary shear wave, sketches may be drawn similar to those presented 
in figure 3, except that the incident wave does not move in the time 
interval 8t. The qualitative character of the downstream disturbance 
is expected to be the same as for problem A. Despite the fact that the 
initial disturbance is not a sound wave, sound waves are produced by the 
interaction and form envelopes when 0 <'l f C 2 or ty cu <1 l f i < n ‘ -*- n 

this case, the values of \jr cZ and t cu are obtained from equation (17) 
with the term involving a^ omitted. The solution, which is presented 
in figure 4, is symmetric about \|r^ = 90°, approaches 90° at M = 1 with 

a half-order singularity, and has the same asymptote at M = « as in 
problem A. As in problem A, when i|r c j <y^<\|f cu , an attenuating pressure 

wave occurs, and for all values of the steady vorticity wave 

appears . 



NACA TN 2879 


13 


The foregoing discussion may be summarized as follows: In case A, 
the incident sound wave refracts, as either a simple sound wave or as a 
more complicated attenuating pressure wave, and an oblique steady wave 
of variable entropy and vorticity appears. In case B, the incident 
sound wave reflects at the shock as a simple sound wave, and a steady 
vorticity wave appears. In case C, the initial "vorticity wave" refracts 
to form a stationary vorticity wave in which, because of the action of 
the shock, the entropy also varies. In addition, a sound wave or pres- 
sure wave is produced by the interaction. 


Solution of Problem A 

In the problem of an oblique plane sound wave overtaken by the 

shock, there are two different solutions to be obtained - one for 

0<ilr, <\|r , and xlr <^,< « and another for i|r ,<\|r..<i|f . The first 

Y 1 T cl cu 1 cl 1 cu 

is the simpler, and will be presented first. 

Solution when 0<i|r-|_-O|r ci or ^ cu <tp< - The initial distur- 
bance ahead of the shock is described in equations (12) and (14). The 
discussion of the preceding section has established that, in these 
ranges of \|r^, the pressure variation behind the shock is associated with 

a simple sound wave. Accordingly, the disturbance pressure is written 

P? (ax-? + Py + a 2 t\ 

i - “ ( r r - -) (18) 

where K, a, 8, and ^ 2 are ‘ to determined. In order that the pres- 
sure satisfy the wave equation (10), 

a 2 + p 2 = 1 (19) 

Equation (18) is written on the assumption that the profile of the pres- 
sure disturbance carries through the shock undistorted, though its ori- 
entation, magnitude, and scale may change. This assumption may be 
regarded as a trial, the correctness of which is inferred from the self- 
consistency of the entire solution so obtained. 

m view of the requirements of the shock relations (equations (5)), 
the arguments of the downstream pressure wave (equation (18)) and the 
initial wave (equation (12)) should match. At the shock, = Vt and 

x„ = (V - U)t. The matching requirement therefore is 

u 
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Equating coefficients of y and t. 



where 

r = U/V (22) 


Equations (19) and (21) yield a quadratic equation for XgAi> the mean- 
ingful solution of which is 



(23) 


where 

n = l/m = tan (24) 

The other solution corresponds to a wave moving in the same direction 
as the shock and is rejected, (inspection of fig. 3(a) shows that two 
families of envelopes might be indicated mathematically, and that only 
one is physically significant.) Thus, the inclination and scale of the 
pressure wave are fully determined, and only the magnitude K remains 
to be found. 

The vorticity wave will also be assumed to have a. profile given by 
the function f. In view of equation (12), its argument can be a func- 
tion only of Xg and y, and must furthermore match the argument of the 

incident wave at the shock, in order that the shock relations may be 
satisfied. At the shock, the argument of the incident wave is given by 
the right-hand side of equation (20). Also, t = (V - U), at the 

shock. The argument of the vorticity wave must therefore be 
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The in cl ina tion and scale of the vorticity wave are thus determined. 


The density fluctuation behind the shock consists of two parts - 
one part associated with the pressure fluctuation to form the sound 
wave, its magnitude determined by equation (ll)j the other part associ- 
ated with the vorticity wave. Thus, 



Likewise, the velocity components each consist of two parts, the 
first associated with the sound wave and the second associated with the 
vorticity wave. Accordingly, 




The requirement that all terms in the shock equations (5) have the 
same functional form and the same argument suggests that 



(Cross-differentiation shows that these two equations are compatible . ) 

The solution is completed by the algebraic determination of the 
various unkn own constants. The coefficients F and H may be found in 
terms of K through equations (7); I maybe found in terms of G from 
the incompressible continuity equation. The remaining unknowns K, Q, 

G, and L may be successively determined by use of the four shock rela- 
tions (equations (5)) when it is recalled that the arguments of all 
quantities have been arranged to match at the shock. Neither the 
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details of the remaining procedure nor of the final solution are partic- 
ularly interesting, and therefore the analysis has been completed in 
appendix B. The numerical results will be discussed in a subsequent 
section. 


Form of attenuating pressure wave when <^i <t cu - - It is 
intended to form the solution for this range of ijr^ in essentially the 

same way as was done in the preceding paragraph. An essential step in 
that solution was the assumption that the refracted sound wave has the 
same profile as the initial disturbance. Therefore, in order to proceed 
with the analysis of the case <tp< t cu it is first necessary to 

determine the form of the pressure disturbance just behind the shock, 
and the manner in which it attenuates with distance behind the shock. 


Tentatively, the pressure is written as a function of two variables 

only: 


where 


Po Pp 

p- = (n>£) (30) 

2 2 


T l= " XJ (^2 • (V - u )t] 

(ccx 2 + Py + cVt) 


(31) 


and the constants d, a, p, and c require determination. The variable 
T) is proportional to distance behind the shock front, -x 2 + (V - U)t. 

Neither the undisturbed shock front nor the incident wave have curvature. 
Therefore it is expected that along any one of a family of planes moving 
with constant velocity, any variation of pressure would be due solely to 
the attenuation associated with distance behind the shock. This consid- 
eration leads to the definition of the second variable £, such that the 
equation £ = constant defines a plane moving obliquely with a constant 
velocity. 


The wave equation 


(10) is satisfied if 



(32a) 
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a 


A a - r)c 
a2^ 


(32b) 


&2 


d 2 (l - r) 2 - c 2 ] = d 2 - (a 2 + p 2 ) 


(32c) 


The boundary conditions to be applied in solving Laplace's equa- 
tion (32a) are 


jr («*£> - 0 (33) 

p 2 

and a condition (equation (5c)) at the shock providing compatibility 
with the initial disturbance. In the following discussion, this infor- 
mation will be used to infer a likely form for the pressure wave. 

Part of the downstream velocity variation is associated with the 
pressure to form an isentropic irrotational flow. For this part of the 
disturbance field, a velocity potential qp(T),0 may therefore be 
defined such that . 


a 


u 2 =c Px 2 = + >^<P£ 


v 2 =c Py =X7 Cp ^ 


(34a) 

(34b) 


In view of equations (7) and (34), 


Jt er * 4 <P - - jr grad p 2 ; g = -<P t 


Therefore, 


gl-Avd'- (34c) 

In the case of the simple sound wave, compatibility at the shock was 
obtained by supposing that the profile of the disturbance (the function 
f) carried through the shock undistorted. In the present case, the cor- 
responding assumption would be that, just behind the shock, at tj = 0, 
the various disturbance quantities each contain a term proportional to 
tit), and that £ should therefore match the argument of the initial 
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disturbance at the shock. From equation (34b), if Vg(0,£) is to con- 
tain a term proportional to f(£), then the potential must contain a 
term proportional to 




T - i 

n 


dt 


which satisfies equation (32a) and hence the wave equation, satisfies 
boundary condition (33), and has the property that (0,£) = f(£).. 

This solution may be regarded as the result of a distribution of singu- 
larities along the plane of the shock (ri =0). In the skew coordinate 
system t|,£, these singularities may be identified as potential- flow 
vortices. Therefore, from equations (34a) and (34c), Ug and pg would 

each contain a term linear in f(£) at the shock and, in addition, a 
term linear in 


poo 

g(O sqp 71 (l) (0,£) = ■- P.V. / (35a) 

where P.V. denotes the Cauchy principal value of an improper integral. 

Of course, V 2 would likely contain a term linear in g(£) at the 
shock also, and therefore the potential would have another part 

f ( t ) In \j \ 2 + (t - t) 2 ~\ dT 


satisfying equations (32a) and (33), and having the properties: 

<^ (2) (o,£) = g(0; <p Tl (2) (o,0 = -f(0 (35b) 

This solution may be regarded as a distribution of potential- flow sources 
along the plane of the shock, in the tj,£ system. 

Thus, the quantities associated with the attenuating pressure wave 
may tentatively be written in the following form: 

^ . k(D (,,{) ♦ *< 2 > h.t) 

r O 



( 36 ) 
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where 








f(T) 


Tl 2 + (T - C ) 5 


dT 


(2) (2) _ m (1) _ 1 

$ = cp^ -cp n - - 


f(T) 


/-co 


T - C 

T ] 2 + ( T- 5) 2 


dT 


(37) 


J 


and the constants k( 1 ) and k( 2 .) require determination. At the shock 
equations (35) and (36) provide that: 

— p— = k(!) f(5) + k( 2) g(U 

F 2 

Examples . - (a) If it happens that f(£) = sin 2n£, then 
g(£) = cos 2 jt£ , and 

= e' 2ltT l sin 2 jc£ 

$ (2) = e" 2,tT l cos 2 jc£ • 

_2 _ e -2jtrj |k( 1) S i n 2n£ + k( 2 ) cos 2ji^J (38) 


Thus, 


and therefore, the disturbance undergoes a phase shift in passing through 
the shock, and subsequently attenuates exponentially with distance behind 
the shock. 

(b) If f(0 = (1 + 6 2 )- 1 , then g(£) = - g(l+ £ 2 ) (see accom- 
panying sketch), and 

% - a , „)* + ? [ k(1) (1 + ^ - k(2) ? ] 
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Solution when i|r„ 7 < i|f.. < \|T .. - The form of the pressure wave has 

-L CU 

been adduced in the previous paragraph (equation (36)). The quantities 
d, a, p, and c may be found by using equations (32) and the requirement 
that the disturbance function f have the same argument ahead of and 
behind the shock: At the shock, x^ = Vt, Xg = (V - U)t, this requirement 

leads to the equations 



Equations (32b) and (39b) yield 



and equation (32c) may be solved for d: 


(39c) 
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The vorticity wave is expected to involve a linear combination of 
the profile functions f and g, just behind the shock, and is not 
expected to change its form subsequently, because it is not time- 
dependent in the x 2 ,y coordinate system and therefore cannot attenu- 
ate, as does the pressure wave, or otherwise change character. A match- 
ing procedure at the shock yields, as previously, the argument given in 
expression (25). 

The remaining analysis parallels that following expression (25) in 
the paragraph Solution when 0 or ^cu < ^l <It ‘ 


_ 1^2 (1) 

R 2 T p 2 


,B -rr£ x 2 - 


m + 


Xl 


+ Q< 2 ) 


g 


l/M 


rr A 2 


xo - ly\ 


\ 


xi 


(40) 


fm + l/M . 

% = F (l) 4 ( 1 ) + P ( 2 ) *( 2 ) + g (D f ( 1 ~ r / 2 + q( 2 ) 

V \ A! 


_1/M 


1 - r 


-2 - l y 


(41) 


^ = h(D 


t 


♦ *<*> + id) f ~ * I< 2 > ~ iy ) 


? - ^ (1) * 


(m + l/M)Vt - l: 

l , t (2) _ 

(m + l/M)Vt - ly 

Xp 


X ! 


(42) 


(43) 


t n / T (l) (m + l/M)Vt - ly ( 2 ) (m + l/M)Vt - ly \ 

‘ h ~ 1 + l/Mm L 1 |_ Xp J S L X x JJ 


(44) 


The various uiiknown constants remain to be determined algebraically 
through equations (7), the incompressible continuity equation, and the 
shock equations (5), as before. These equations suffice to determine a 
greater number of constants than were required in the previous case 
because the functions f, g, '$(!), and <|>( 2 ) and their derivatives form 
two separate groups of functions whose coefficients may be separately 
equated. Details of this procedure are provided in appendix B. 
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Solution of Problem B 

The analysis of problem B (sound wave overtaking shock from behind) 
is identical in all essential respects to that of problem A when only a 
simple refracted sound wave is involved. The only differences which 
arise are the slightly different matching of arguments at the shock, and 
a slightly different form assumed by the shock relations (5). 

The equation (18) is adopted in the present case to represent the 
reflected sound wave. Matching the argument of this expression with that 
of the initial disturbance (equation (14)) yields 



Equations (19) and (46) provide a quadratic equation for ^22^21* the 
useful solution being: 


*22 

*21 



(46c) 


The other solution is X 22 A 2 I = corresponding to the incident wave 
itself. 

The right side of equation (45) is used as the argument in the 
expressions for and | y.: 
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The right side of equation (45) also yields the argument of the vorticity 
wave, when the substitution t = Xg/ (V - U) is made: 


a 2 1 


L m + vTT" ?)J 

x 2 - 2 y 


^21 


(48) 


Equations (26), (27), and (28) maybe adopted to complete the description 
of the flow, except that expression (48) must be used for the argument 
of the vorticity wave. The analysis is completed in appendix B. 


Solution of Problem C 

The only differences between problem A and problem C (stationary 
vorticity wave overtaken by shock) involve the shock equations (5), and 
the matching of arguments at the shock. The difference in matching is 
due to the fact that in the present case, the disturbance is stationary, 
while in problem A, the disturbance moves with velocity a^. Accord- 
ingly, when 0 < tyq < equations (18), (21), (23), and (25) through 

(29) apply directly to the present case, if the quantity l/M is omitted 
wherever it appears explicity. When \|r 7 < i|r.. < i|r , equations (36), 

C L _L CU ** 

(37), and (39) through (44) may also be adopted, again omitting terms 
proportional to l/M. The remaining details of the analysis are provided 
in appendix B. 

As previously mentioned, an initial disturbance of type C may have 
a third fluctuating velocity component parallel to the shock, which might 
be represented as follows: 



(see equations (15)). The amplitude A 5 is arbitrary, within the limi- 
tations of linear analysis. This disturbance will pass through the shock 
unaffected, and become part of the steady vorticity wave behind the shock. 
Thus, 
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The shock front itself will not be affected by this purely transverse 
disturbance. 


Singularities at i|r . and y 


Equation (23), whiqh applies in the refracted sound wave solution, 
contains a radical which vanishes at \|» c ^ and \|r cu , and becomes imaginary 

when < < ty cu - ^ be shown that the quantity under the radi- 

cal vanishes with a nonzero slope. Accordingly, the quantity ^2/Al, 
though finite, has a half -order singularity at \| r c j and i|r cu . Further- 
more, this quantity is involved in all the formulas characterizing the 
refracted sound wave (see appendix B). In the range T|r cl < tq < t C u> 

the attenuation coefficient d (equation (39d)) vanishes with half- 
order singularities at if cl arid \|r cu , and similarly affects the remainder 

of the analysis. 


The reason for this singular behavior may be inferred from figure 3. 
According to figure 3(a), when 0 < \|/q < if c j, the refracted sound wave 


is the envelope of an imagined succession of cylindrical waves, as shown 
in the sketch. At \|r c j, however, according to figure 3(b), the cylindri- 
cal waves all meet at a common point of tangency, as shown in the sketch. 
Therefore, successive waves reinforce at one point, giving a singularity 
of the flow. This singularity of course depends on the fact that the 
theory is linear. An exact analysis would presumably show steep though 



This situation is similar to that arising in the linearized analysis 
of compressible flow about bodies: as the Mach number of the flow 

approaches 1, the Mach waves have a common point of tangency at the nose, 
and the wave drag shows a reciprocal half -order singularity in Mach num- 
ber. In the present case, the physical quantities remain finite, but 
have infinite rates of variation with 
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The preceding discussion applies to problems A and C, but not to B. 


RESULTS AND DISCUSSION 

In the following paragraphs, the results (presented in graphical 
form) will be described for problems A and B. Results for problem C are 
presented in reference 6 . The solution for each problem has essentially 
three elements : 

1. Disturbance of the shape of the shock front 

2. Characteristics of the isentropic pressure wave behind the shock 

3. Characteristics of the steady vorticity wave behind the shock. 

Computations have been carried out for three Mach numbers, 1, 1.5, 
and infinity. Of course, the case M = 1 is. really degenerate, because 
the shock is then a weak sound wave, and therefore the interaction with 
the incident disturbance is obtained by linear superposition, the initial 
disturbance passing through the "shock" with no change. Also, M = 1 is 
a singular point, because the range of angles of incidence of initial 
disturbance for which the attenuating pressure wave appears in problems A 
and C vanishes with a half -order singularity (fig. 4), and in problem B, 
because the incident wave is unable to overtake the "shock" moving with 
sonic velocity. 

The results will show that the critical angles .and \[r cu are 

also singular points. In many instances there are not a sufficient number 
of points near the singular points for which computations have been made, 
so that not all of the curves can be faired with complete confidence. 

For this reason, the computed points are shown as heavy dots, so that the 
basis for the fairing will be clear in each case. 

No mention will be made of the temperature disturbance behind the 
shock, which may be obtained directly from equation (4) if the pressure 
and density disturbances are known. 

All disturbance quantities found by the linearized analysis of the 
present report will be proportional to the intensity of the incident 
wave. Therefore, results are divided by the pressure amplitude of the 
incident wave, A 3 (see equations (l2a) and (14a)). . 


Problem A - Shock Overtaking Sound Wave 


1. Shock front disturbance . - In figure 5 are shown the amplitudes 
1 = lS^/A$ an^l"~" 2 )/. A 3 of the incremental velocity of the 
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shock front, due to the interaction. From equations (29), (43), and (44) 
these amplitudes are associated, respectively, with the functions f 
(which defines the profile shape of the incident wave) and g (which is 
an additional profile function arising when t c Z < ♦£ <t C u) to § ive the 

actual incremental velocity. The variations with ^ are quite extreme, 

particularly at the critical angles, where, in fact, there are half -order 
singularities. The variations with Mach number are equally severe. 

When \|ri = 0 (incident wave moving parallel to and toward the shock) 
the figure shows that the shock front is retarded by a pressure wave. 

In the case M = 1, this is because the velocity in an incident compres- 
sive sound wave, relative to which the shock (really a weak compression 
wave) propagates as it moves through the disturbance, is directed against 
the shock. When the incident compression wave moves in the same direc- 
tion as the shock (\|r^ = jt), the shock front is speeded up for low shock 

Mach numbers; at M = 1, this is true because the incremental velocity 
due to the incident wave is in the same direction as the ” shock” movement. 
Whether i|rq = 0 or jt, there is a smaller accelerating effect due to the 
higher velocity of sound in the incident compression wave - that is why 
the curve for M = 1 is not perfectly anti-symmetrical about = Jt/2. 

When M = », the curves are symmetrical about = tc/ 2 because 

the incremental sound wave velocity is vanishingly small compared with 
the shock velocity. 

For each value of M there is a value of ^ for which the shock 

intersects the incident sound wave permanently at one point on the 
traveling wave, and the problem becomes essentially steady, so that the 
increment in shock velocity vanishes (though a steady displacement 

occurs). This happens when a^ = - mV, or, \|f]_ = sec“^(-M), yielding 

131.8° when M = 1.5, and 90° when M = oo, and is the case of steady 
interaction of a Mach wave and a normal shock which has been treated by 
Adams (reference 2). 

2. Characteristics of pressure wave . - (a) Inclination of refracted 
pressure front: Figure 6 shows the angle between the directions of 

propagation of the shock and the refracted pressure wave behind the 
shock. In view of equations (18) and (31), this quantity is given by the 
equation 


if 2 = cot’ 1 (“ a /P) 

Of course, when ^ c z < if± < this wave is not a sound wave, and the 

inclination shown refers to a front parallel to which physical quantities 
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depend only on distance behind the shock. Outside the range 

^cl < tl < tcu> the pressure wave is a sound wave. When M = », the 

curve is symmetrical about tyq = «/2. 

(b) Coefficient of attenuation of pressure wave: Figure 7 shows 

the coefficient d which appears in equations (31) and (39d). This 
quantity vanishes at \}r c ^ and \|r cu , indicating that the sound wave 

solution and attenuating wave solution meet continuously at the critical 
angles. When M is either 1.5 or the maximum value of d is 

about 1. This implies a rather rapid attenuation - if f(£) = sin 2jt£, 
it has been shown (equations (38)) that the wave attenuates as 
exp( -2jtT)) . From the definition of q (equation (31)), when d « 1, the 

attenuation factor becomes e - '*' at a distance behind the shock approxi- 
mately equal to l/2jt times the wave length of the incident sound wave. 

(c) Propagation velocity of pressure wave when \|r C 2 < \|rq < ty cu : 

Figure 8 shows the quantity c of equations (31) and (39c), combined 
with other quantities to give propagation velocity as a fraction of the 
speed of sound ag. Since the solution when < ^1 < ^cu mee ts the 

sound wave solution continuously, the propagation velocity is equal to 
the speed of sound at \|; c ^ and \|r cu . The change of sign of c is 

taken into account in figure 6 by the 180° shift of direction shown at 
the angle for which c = 0. 

(d) Ratio of scales of pressure waves behind and ahead of shock: 
Figure 9 shows the quantity Ag/Xq (equation (23)), which was defined 

only for the sound wave solution. However, inspection of equation (31) 
shows that the equivalent quantity when t|r c j < tyq < \[r cu is 

(a 2 + p 2 ) -1 / 2 . 

The reversal of sign of Xg/Xq signifies a reversal of the direc- 
tion of propagation of the pressure wave relative to the shape of the 
incident wave. For example, when \|r^ = 0, the incoming and outgoing 

waves might appear as follows: 




/V 


incident wave profile 
Xg/Xq positive 


Shock 
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whereas when i|r^ = it, they would appear thus: 


negative 

The difference between these two cases consists of a difference in sign 
in the arguments of the refracted wave in the two cases, and arises 
formally in the present analysis as a change in sign of X 2/^1* 

When M = 1.5, th$ magnitude of Ag/^l is S rea -ter for ^ = it 
than for ^ = 0, because when the incident wave is traveling in the 
same direction as the shock, the shock requires a longer time to traverse 
the incident wave than when the two waves travel in opposing directions. 

(e) Amplitudes of pressure disturbance behind shock: Figure 10 

shows the coefficients (of f, or of $(D and $(2) when 

^cZ < +1 < "which describe the pressure wave behind the shock; 

k(!> 2 ) s (see equations (18) and (36)). As in figure 5, the 

flow is shown to be singular at the two critical angles and to vary 
markedly with both ^ and M. At M = 1.5, near ^ = it, the 

refracted sound wave is seen to be very weak. . 

Since the pressure wave is isentropic, the coefficients of the cor- 
responding part of the density variation are equal to l/y* times the 
pressure coefficient. 

(f) and '(g) Coefficients for velocity components in pressure wave: 
Figures 11 and 12 show the coefficients for the velocity components 
associated with the isentropic pressure wave behind the shock (equa- 
tions (27), (28), (41), and (42)). The velocity resultant is longitu- 
dinal with respect to the direction of propagation of the pressure 
wave, except when ] < \|fq < ^ C u* 

3. Characteristics of the steady vorticity wave behind the shock: 

(a) Inclination of the steady vorticity wave: Figure 13 shows the 

inclination ^3 of the vorticity front behind the shock. From equa- 
tion (25), 

= cot " 1 




Shock 


NACA TN 2879 


29 


(b ) Ratio of scale of vorticity wave to that of incident sound 

wave: From equations (12a) and (25) this is 

(i 

which is shown in figure 14. When M = 1 and tyq = Jt , the scale ratio 
goes to 00 because the "shock" is unable to overtake the incident wave 

(c) Amplitude of density variation in vorticity wave behind shock: 

Figure 15 shows ' and ' (equations (26) and (40).). Appar- 

ently, when M is of order 1.5, the vorticity wave is very weak. 

(d) and (e) Coefficients for velocity components in vorticity wave 
Figures 16 and 17 show the coefficients for the transverse velocity 
field associated with the vorticity wave. 


Problem B - Sound Wave Overtaking Shock from Behind 

Consistent with the previous analysis , computations have^been 
carried out only for values of 4^1 (incident wave inclination) suffi- 
ciently close to 180° so that the component of propagation velocity in 
the direction of the motion of the 'Shock is greater than the velocity 
of the shock relative to the fluid behind. These values of typ are 

134.5° and 112.2° for M = 1.5 and co, respectively. 

1. Shock-front disturbance . - Figure 18 shows the amplitude of the 
incremental shock-front velocity, L' = L/A 3 . When 1^21 is near. 180°, 

an incident pressure wave displaces, the shock ahead. 

2. Characteristics of reflected sound wave . - The downstream pres- 
sure wave in problem B is always a simple sound wave. 

(a) Inclination of reflected wave: The inclination 

^22 = cot"-*-(-a/p) (see equation (18)) is shown in figure 19. In effect 

the incident and reflected waves coalesce into a single wave at the 
critical angles. 

(b) Ratio of scales of reflected and incident sound waves: This 

ratio is given by equation (46c) and is shown in figure 20. At \|rgi = 
the scale ratio is greater when M = 1.5 than when M =», just as in 
problem A, and for the same reason. At the critical angles, the ratio 
becomes 1 , because the incident and reflected waves coalesce. 
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(c) Amplitude of reflected sound wave: The pressure anplitude 

K' = K/A3 (equation. (18)) is shown in figure 21. At = it, an inci- 

dent compression wave reflects as an expansion wave with a strength 
which is greater the higher the Mach number, but always less than that 
of the incident wave. 

The velocity components associated with the reflected sound wave 
are obtained simply from K' and are therefore not plotted. 

5. Characteristics of the steady vorticity wave . - The coefficient 
Q 1 of the density fluctuation in the vorticity wave is simply propor- 
tional to L' (combining equations (B26) and (B27) of appendix B) and 


is therefore not shown in 

a figure 



At 

M = 1.5; 

Q' 

= - 0.147L ' 

At 

3S 

M 

8 

Q' 

= - 1.43L' 


(a) Inclination of the vorticity wave : From equation (47), the 

inclination is given by 


+23 ’ cot ' 1 ( 5 1 + 


a 2 /V 

m(l - r) 


which is plotted in figure 22. 

(b) Ratio of scale of vorticity wave to that of incident sound 
wave: As in problem A, this ratio is given by (l /l) sin and is 

shown in figure 23. As for the reflected sound wave, the scale ratio 
is larger for M nearer one. 


(c) and (d) Velocity' variations in vorticity wave : Coefficients 

G' and I' of the transverse velocity fluctuation (equations (27) and 
(28)) in the vorticity wave are shown in figures 24 and 25. 


4. Wave reflection at the critical angle . - The analysis and fig- 
ures show that at the critical angle the incident and reflected sound 
waves coalesce to form a single sound wave. This statement may be inter- 
preted to mean that a sound wave incident at the critical angle reflects 
as a steady vorticity wave only. Therefore, the shock disturbance and 
vorticity wave characteristics may be expressed in terms of the pressure 
anplitude . of a single incident sound wave of strength 
A 3 + K = A 3 (l + K’). 
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CONCLUDING REMARKS 

In principle, the interaction of a shock with any weak flow field 
may he obtained by first constructing the initial flow field as a 
linear combination of plane waves of varying strength and orientation. 
From the present analysis the interaction of each constituent wave with 
the shock may be found. Assembling the resulting waves behind the shock 
then would yield the desired solution. Regrettably, the formulas for 
the interaction depend on the angle of incidence in a rather complicated 
way and it would in general be difficult to evaluate explicitly inte- 
grals in which these formulas are used for the distribution functions. 
Numerical procedures could be used for this purpose, though a technique 
would be required for dealing with the singularities at the critical 
angles \|r cl and \lr cu . 

The nature of the solution of problem B perhaps requires clarifi- 
cation, in that the angles of incidence have been restricted to 

the range for which the sound wave overtakes the shock. For purposes 
of superposition, however, all incidence angles must be considered. 
Inspection of figures 6 to 12 shows that a sound wave of any incidence 
angles between 0 and it may be identified either as an incident 
wave or as a reflected wave, in the sense of the analysis. For purposes 
of linear superposition, the distinction between incident and reflected 
waves is of no . significance. The point of view may be adopted that 
whatever its angle, a constituent wave of the initial flow has asso- 
ciated with it another sound wave of another angle and a steady vor- 
ticity wave. The notion of cause and effect is not needed. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, October 15, 1952 
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APPENDIX A 
NOTATION 

The following symbols are used in this report: . . . 

Al,A 2 ,A 3 ,A 4 coefficients of incident disturbance (equations (12a), 
(14a), or (15)) 

a velocity of sound 

coefficients in shock relations (equations (6)) 

C^_,C 2 coefficients in shock relations (equations (6)) 


Dl,D2 

d 

F 

f 

G 

g 

H 

*n 


dimensionless velocity of propagation of pressure wave 
behind shock (equation (31) or (39c)) 

specific heat at constant volume 

coefficients in shock relations (equations (6)) 

coefficient of attenuation of pressure wave behind shock 
(equation (31) or (39d)) 

coefficient of part of u 2 associated with pressure wave 
(equation (27) or (41)) 

profile function of incident wave (equation (12a), (14a), 
or (15)) 

coefficient of part of u 2 associated with vorticity wave 
(equations (27) or (41)) 

additional profile function appearing behind shock 
“ P.V. f (t) (t - £)”■*■ dT , equation (35a) 

coefficient of part of v 2 associated with pressure wave 
(equations (28) or (42)) 

functions involved in solution when i| < ^ < i|r 
(equations (B19), appendix B) 
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I 

J 

K 

L 

Z 

M 

m 

n 

P 

P 

Q 

R 

r 

t 

U 

u 

V 


v 


x l> x 2 


y 


coefficient of part of V 2 associated with vorticity wave 
(equations (28) or (42)) 

gas constant (equation (4) ) 

coefficient of pg (equations (18) or (36)) 

coefficient for £(y,t) (equations (29) or (43)) 

sin \|t-l or sin 

Mach number of shock (= V/a-^) 

cos or cos \|f2i 

tan or tan ijr g ^ 

mean static pressure 
perturbation in static pressure 

coefficient for density fluctuation in vorticity wave 
(equations (26) or (40)) 

mean gas density 

ratio of U to V 

time 

mean velocity of gas behind shock (fig. l) 

perturbation of velocity component in x-direction (fig. l) 

mean velocity of propagation of shock in gas at rest 

(fig. 1) 

perturbation of velocity component in y-direction 

coordinates measured in the direction of the shock propa- 
gation, relative to which the gas is (on the average) at 
rest ahead of and behind the shock, respectively 
(figs. 1 and 2) 

coordinate orthogonal to xq or X2 (figs. 1 and 2) 
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<x,p 

r 

i 

n 

0 

8 

X 

t 

P 

°,x 

$ 

<T> 

tl^2i 

^ 2^22 

+ 3 '+23 


functions defining pressure front (equations (21) or (39)) 

ratio of specific heats (*1.4 for air) 

variable upon which pressure wave depends when 
tel < +i < ^cu (equations (31)) 

variable upon which pressure wave depends when 
V cl < t 1 < i cu (equations (31))' 

mean static temperature 
fluctuation, in static temperature 
scale of pressure wave (equation (23)) 
displacement of shock front (fig. l) 
fluctuation in gas density 

functions appearing in solution for L (equations (B9) or 
(B3l), appendix B) 

function associated with pressure wave when 
^ cl < t 1 < t cu (equations (37)) 


velocity potential associated with pressure wave when 

\|r , < < \|r 

cZ Y 1 Y cu 

angles of inclination of incident waves (fig. 2) 

angles of inclination of pressure wave behind shock (figs. 6 
or 19) 


angles of inclination of vorticity wave behind shock (figs. 13 
or 22) 


^cZ’^cu lower and upper bounds, respectively, of the range of 

for which the attenuating pressure wave appears (equa- 
tion (17)) 

Subscripts: 

Subscript notation for partial differentiation has been used where 
convenient . 
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1 conditions ahead of shock 

2 conditions behind shock 

Double subscripts 21 and 22 (problem B), incident and reflected 
waves, respectively. 

Superscripts : 

(1) coefficients associated with f at' the shock 

( 2 ) coefficients associated with g at the shock 

Primed coefficients are referred to intensity of incident wave. 
(A 3 in problems A and B, in problem C). 

Primes are also used to denote ordinary differentiation of f 
and g with respect to their arguments. 
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APPENDIX B 


COMPLETION OF INTERACTION ANALYSIS 
Problem A 


The solution of problem A is found in two separate ranges of 

for which the solution involves a refracted sound wave in one case and 
a refracted attenuating pressure wave in the other. The analysis will 
be completed in that order. 

1. 0 < \(r < \|r , , ^ <+_<«.- Equations (27) and (28) are sub- 

^ Co CU 

stituted into equations (7), and coefficients of f' are equated, 
yielding : 


> (Bl) 


When equations (26), (27), and (28) are substituted into continuity 
equation (8), the quantities associated with the sound wave combine to 
satisfy continuity. The terms associated with the vorticity wave then 
must satisfy the incompressible continuity equation (i^ + Vy = 0) 

because the corresponding density term is time -independent. Therefore, 
equating coefficients of f ' , 


F = 


a 2 

- Ka 
yV 


H = 

V 


rv 


Kp 


1 + l/Mm _ 

= n(i Try G 


(B2) 


The shock equations remain to be satisfied; it should be noted that the 
arguments of f for the various quantities have been matched at the 
shock, so that coefficients of f maybe equated. From equations (5a), 
(27), (29), and (12a), 

G = L - F - B-l(L - Aj) - B 2 A 3 (B3) 


From equations (5b), (5c), (5d), respectively, with the necessary sub- 
stitutions made, 


Q = C 1 (L - A-l) + C 2 A 3 - ^ 


(B4) 


K = Di(L - A x ) + D 2 A 3 


(B5) 
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L = 1 ± y m (H + I - k 9 ) ' (B6) 
nr • c 

The terms H and I in equation (B6) may he expressed in terms of L 
by using equations (Bl), (B2), (B3), and (B5). The resulting equation 
for L yields the result 



or, from equations (6), (21), and (23): 



and where use is made of the results from equations (2), (6) and (22) 
that 


r + i 




(i 


r) 



r 
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when L' is known, the remaining coefficients follow from equations (Bl) 
through (B5) 

2. ty cl < ^ < i|f cu . - Equations (36) and (4l) substituted into (7) 

give 

F (i) $^_(l) + f (2) «j^( 2) _ _ _JL_ Af(i) <^(l) + $ x ( 2 )^ (BIO) 

From equations (31) and (37), 


„ (1,2) = (V - U) ** (1,2) +$V $ (1,2) 
x a x n Ai s 


S (1 ' 2) "^V 1 ' 2)+ ^V 1 ' 2) 

Equation (37) show that 


S> 


(1) _ * (2) (1) _ . (2) 

1 ' #£ -’®T1 


Accordingly, equation (BIO) may be written: 

(1 - r) dP(l) * (1) - cf(D $ (2) + (1 - r) dp(2) $ (2) + cF (2) . (l) 

T\ T] ^ T) » T) 

T dK ^ - dK^ 2 ) fc^ 2 ) + aK^ 2 ) $ TJ ( 1 )'| 

Equating coefficients of and (J^ 2 ) separately, 


(1) and 

$ 

n 

n 

+ cF^ 2 ) = 

a 2 2 


rv 2 


a 2 

1 

O 

II 

2 


r v 2 


(Bll) 
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Similarly, equations (36), (42), and (7) yield: 


(l - r) cIh(1) + ch(2) = - pK( 2 ) 


(1 - r) aH( 2 ) - ch(D = Pk(D 

rv 2 


When equations (41) and (42) are substituted into equation (8), and 
coefficients of f } and g 1 are equated, the following equations are 
obtained; 


Mm G (l) 


nil - r, 


1_ 

__JSL(}(2) 

n(l - r) 


The shock conditions remain to be applied. Equations '(5a), (5b), 
(5c), and (5d) yield, respectively, the following four pairs of equa- 
tions , when the coefficients of f and g are separately equated: 


G^ 1 ) = i/ 1 ) - - B^ (i/ 1 ) - Aj] - B 2 A 3 


g (2) _ l (2) _ f ( 2) _ Bl x,(2) 


q^ 1 ) = cmi/ 1 -) - a-J + c 2 a 3 - K^Vr 


Q( 2 ) = Cnl/ 2 ) _ K^ 2 Vr 


k( x ) = d-l (b 1 ) - a-A + d 2 a 3 


k(2) = D]L( 2 ) 
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L (l) _ 1 + l/Mm ^j(i) + j(l) _ 

L (2) = 1 ± / h ( 2) + x(2)j 


(B17) 


Equations (Bll) through (B17) may he combined to yield the solu- 
tions for l(D and l( 2 ) 


L 


(D’ = 



^3 + hj^ 


L(2)' = 



(B18) 


where 



h-z = - 


2d(l - r) 


ia ( x + HH 1 r ) 




(B19) 


Equations (B16) are substituted into equations (Bll ) , which may be solved 
to yield 


F W -_ h4 ( L U)' + i + |). h 5 L( 2 )' 

F (2>' =h 4 L(2>' + h 5 ( L a)’ + J + g 


(B20) 
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where 


* Y + 1 (1 + o 2 ) 

_ 2rdo(l - r) V 2 
7(1 + a 2 ) ■ ap 2 


Similarly, 


= (1) '=^4 (L(l)' + iL + 5 5) + o h5 L( 


h( 2 ) ’ = ~ h 4 L^ 2 ^ 
a * 


( L d)' + jl + 

V rM d-l 


The remaining quantities follow from equations (B14) and (B15), 


Problem B 

The solution of problem B is restricted to angles for which 

V , . 

-m > — (1 - r ) , and involves a reflected sound wave and vorticity wave 
a 2 

behind the shock. The momentum equations (7) provide 


- — oK 

rv 


a 2 

- -I PK 


and the continuity equation (8) provides that 


_ i / a 2 / V 

I = — m + - 

l V 1 - r 
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All perturbation quantities vanish ahead of the shock, and shock rela 
tions (5a), (5b), (5c), and (5d) yield, respectively, 


G = L(l - B x ) - F - A-l 

(B25) 

Q = CjL - a 4 - K/r 

(B26) 

K = D-jL - A 3 

(B27) 

m(l - r) + a ? /V 

vr^- (H + 1 + a 2 } 

(B28) 


Solving the foregoing set of equations for L, 


, = _l = r -n L + 2_z_i r \ 2L_ _ 

A3 * ' 2 ' l + X-!^r(l +r 2i:) 


(B29) 


where 




O =s 


1 +■ 


n 

-^7v" 


m 


(1 - r) 


= - ( Zcr + m) 


1 - r 


1 + ]L 4- i r y 


(B30) 


The remaining quantities follow from equations (B23) through (B27) 


Problem C 

The analysis of this problem parallels that of problem A. 

1. 0 < ^2. < i|r cJ , t cu < tq < - The equations of motion (7) and (8) 

provide 



F = 
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The shock equations (5) give 

G = L - F - B-^L - A ± ) 

Q = C 1 (L - A x ) - K/r 

K = D 1 (L - A ± ) 

L = i (H + I - A 2 ) 


These equations may be solved to yield 



2. < ^ cu * - The analysis is essentially identical to the 

corresponding part of problem A, except that the constants B^, Cg, and 
Dg do not appear in the present case, and the Mach number M does not 

appear explicitly. Except for these differences, equations (B14) 
through (B17) may be carried over to the present case. The solutions 

for and may be written as 
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L (l)' s L (l) h i h 2 + Zh 5 2 

^AjZ + A 2 2 h 2 2 + h 3 2 

l ( 2) ' s L^ 2 ) = h3 h l ~ Zh 2 

Va-l 2 + A 2 2 h 2 2 + h 3 2 


The definitions of the h's follow those of problem A (equations (B19) 
and (B2l)), with the exceptions that 


h 


1 = 



and, in this case, a = n(l - r). The solutions for the F's and H's 
are : 



f ( 2 )’= l (Z) 

V A 1 2 + A 2 2 

H (l)’= h(1) = 

V A i 2 + A 2 2 

V A 1 2 + A 2 2 


= h 4 L^ 2 ) ' + h 5 (i/ 1 ) ' - i) 

i h 4 (l(D ’ - z) + cxh 5 L( 2 )' 

i h 4 L( 2 ) ' - crh 5 (i/l) ' - z) 


The remaining quantities follow as before. 
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no disturbance 

Figure 1. - Notation for shock wave propagating into a region of weak 

disturbance . 



(a) Shock overtaking sound wave. (b) Sound wave overtaking shock from behind. 



Figure 2. - Types of initial disturbance considered. 
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Figure 3. - Formation of waves behind shock, because of interaction with 

sound wave (problem A . ) 



Critical angle 
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Shock Mach number, v/aj_ = M 

Figure 4. - Critical angle of incident wave for formation of refracted sound 

wave behind shock. 







Propagation velocity of pressure wave, — . JL Attenuation coefficient 

_y a 2+p2 & 2 
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Figure 7. - Problem A:, attenuation of pressure wave behind shock (shock 

overtaking sound wave). 



60 70 80 90 100 110 120 130 140 150 160 170 180 


Inclination of incident wave, deg 

Figure 8. - Problem A: propagation velocity bf pressure wave behind shock (shock 

overtaking sound wave). 




Ratio of pressure wave scales, 
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Figure 9. - Problem A: ratio of scales of pressure wave up- and down- 

stream of shock (shock overtaking sound wave ) . 




Coefficient of u ? associated with pressure wave 





Coefficient of vg associated with pressure wave H 
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Inclination of incident wave, deg 


Figure 12. - Problem A: coefficient for part of v 2 associated 

with pressure wave (shock overtaking sound wave.) 




Inclination of vorticity wave, ^ 3 , deg 
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0 20 40 60 80 100 120 140 160 180 


Inclination of incident wave, \h, deg 

Figure 13. - Problem A: Inclination of vorticity wave behind shock 

(shock overtaking sound wave.) 



Ratio of scales of vorticity wave and 
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0 20 40 60 80 100 120 140 160 180 

Inclination of incident wave, tyg, deg 


Figure 14. - Problem A: ratio of scales of vorticity wave behind 

shock to that of incident sound wave (shock overtaking sound 
wave . ) 







Coefficient of uo associated with vorticity wave , G 
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Figure 16. - Problem A: coeff 

with vorticity wave (s' 



Coefficient of Vg associated with vorticity wave. 








Inclination of reflected sound wave, ^ 22 ’ ^- e .6 
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I ‘ 1 1 i 1 , 1 1 1 I 

90 100 110 120 130 140 150 160 170 180 

Inclination of incident wave, ^ 21 * deg 


Figure 19. - Problem B: inclination of reflected sound wave 

(sound wave overtaking shock.) 
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Figure 20. - Problem B: ratio of scales of reflected and incident 

sound waves {sound wave overtaking shock.) 




Pressure amplitude of reflected sound wave 
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Figure 21. - Problem B: amplitude of reflected sound wave 

(sound wave overtaking shock.) 


i 
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Figure 22. - Problem B: inclination of vorticity wave produced 

behind shock (sound wave overtaking shock. ) 



Inclination, of incident wave, tgi* deg 


Figure 23. - Problem B: ratio of scales of vorticity wave behind 

shock and incident sound wave (sound wave overtaking shock.) 



66 NACA TN 2879 



with vorticity wave (sound wave overtaking shock.) 



Figure 25. - Problem B: coefficient for part of v 2 £ associated 

with vorticity wave {sound wave overtaking shock.) 
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